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I. INTRODUCTION 

Currently, heavy-to-heavy form factors are mainly derived from models. Since models 
are used in Monte Carlos, which in turn feeds into studies of backgrounds and efficiencies 
in many experimental settings, constraining them is very important. Model independent 
information on the form factors can be obtained from the Heavy Quark Effective Theory 
(HQET). HQET [I|-0 vastly simplifies the nonperturbative calculation of form factors by 
relating all of them in the infinite mass limit to a single universal Isgur-Wise function, ^, 
which describes the nonperturbative physics of the light degrees of freedom in the heavy 
mesons. This function is normalized to unity at zero recoil, where the heavy meson in the 
final state has the same velocity as the initial one. Nonetheless, HQET does not predict the 
explicit form of the Isgur-Wise function. 

Since it is not possible to calculate heavy-to-heavy form factors from first principles, 
the next best thing to do is to theoretically bound them.Q A set of model independent 
bounds on form factors have been derived |^,§|, and can serve as a consistency condition for 
phenomenological models. Bounds to 0{l/mQ) at arbitrary momentum transfer had been 
presented in 0J^ for heavy-to-heavy and heavy-to-light form factors. However, as noted in 
0-^, the lower bound is sensitive to the perturbative corrections, and leading order in a^ 
correction should be included to provide more reliable bounds. Therefore, in this work we 
present an analysis of the leading QCD corrections on the bounds of individual heavy-to- 
heavy, more specifically B -^ D^*^iu, form factors. The information of these bounds will 
help us rule out unrealistic models, and may in turn be applied to certain decay amplitudes 
that are of interest. 

This paper is organized as follows. In Section II, we review the derivation of the sum 
rules used to obtain the model-independent bounds on form factors defined in Section HI. 
Section HI also gives the proper combinations of structure functions used later for the 



^See, however, 0] for model independent paranieterizations of the form factors. 



bounds on each of the form factors. In Section IV, we perform the tree level nonperturbative 
expansion and first order QCD perturbative calculation to eliminate the leading uncertainty 
in the bounds. Section V provides the bounds on individual form factors explicitly with the 
structure functions given in the Appendix, and discusses various influences of the parameters 
appearing in the expansion of the bounds. Some popular form factor models are compared 
with our bounds in Section VI. Order a|/5o corrections to the bounds at zero recoil are 
computed in Section VII. The conclusion of our work is summarized in Section VIII. The 
Appendix lists the perturbative corrections to the structure functions. 

II. REVIEW OF SUM RULES 

The sum rules are derived by relating the inclusive decay rate, calculated using the 
operator product expansion (OPE) and perturbative QCD, to the sum of exclusive decay 
rates. In the following, we follow P,H in the derivation of the bounds. 

First, consider the time ordered product of two currents between B mesons in the mo- 
mentum space, 

T'' = ^J d'^ ^"'"^ iB{v) I TiJ^\x)rm \Biv)) 

= -g^^T, + v^^v^T^ + le^^-^q^VpT, + q>^q^T, + (g^t;^ + v>^q'')T,, (1) 

where J'^ is a 6 — *> c axial or vector current. The time ordered product can be expressed as 
a sum over hadronic or partonic intermediate states. The sum over hadronic states includes 
the matrix element {H\J\B)^ where H = D,D^*' is the final state heavy meson of interest. 
After inserting a complete set of states and dotting with four vectors a*a,y, one obtains 

+^EP.)»*»fe-,-)Jf^^g^. (2) 

where T(e) = a*^T'^^ai,, e = Mb — Eh — v ■ q, and the sum over X includes the usual 
d^p/2Ex for each particle in the state X. We choose to work in the rest frame of the B 



meson, p = Mbv, and the z axis pointing in the direction of q. We will hold qs fixed while 



analytically continuing t> ■ g to the complex plane. Eh = y Mfj + gf is the H meson energy. 
There are two cuts in the complex e plane, < e < oo, corresponding to the decay process 
b ^ c, and — oo < e < —2Eh, corresponding to two b quarks and a c quark in the final 
state. The second cut will not be important for our discussion. 

The integral over e of the time ordered product, T{e), times a weight function, e"WA(e) 
can be computed perturbatively in QCD P,|l. For simplicity, we pick the weight function 
Wa{^) = ^(A — e), which corresponds to summing over all hadronic resonances up to the 
excitation energy A with equal weight. Relating the integral with the hard cutoff to the 
exclusive states requires local duality at the scale A. A must be chosen large enough so that 
the structure functions can be calculated perturbatively. 

Taking the zeroth moment of T(e) 

2m Jc ^ 2Mb 

^ \{H[v')\a.J\B{v))f 

AMbEh ^ ' 

gives an upper bound on the matrix element. 
The first moment of T(e) gives 

^.fdeee{A~e)T{e) 

= Y.O{A-E^ + En) {E^ - E^) {2.f6\p^ + g) '^^' "/ J^^'' 

X^H AMbEx 

> {E, - En) E d{A -Ex + En){2n)^'{px + g) '^^' ?/ J^^'' - (4) 

X^H AMbEx 

where Ei denotes the first excited state that is more massive than H meson. Here the 
validity of the inequality relies on the assumption that multiparticle final states with less 
energy than Ei contribute negligibly. This assumption is true in large Nc-, and is also 
confirmed by current experimental data. 

A lower bound can be formed by combining Eqs. (§) and @) to be 



iH,n\a-j\Bm\^^^,^^_^^,^^^^^ (,, 



AMbEh -2mJc " ' " ' \ El -Eh 

Therefore, we find the bounds 

2711 Jc ^ ' ^'\ El- Eh) ~ AMbEh 

As emphasized in [^, the upper bound is essentially model independent while the lower 
bound relies on the assumption about the final state spectrum. These bounds can be used 
for the decays at arbitrary momentum transfer g^ and are good for both heavy mesons and 
baryons. (For baryons, a spin sum ^^ J2s,s' needs to be included in front of the bounded 
factor.) 

Since l/{Ei — Eh) ~ I/Aqcd, the lower bounds will be limited to order in I/tjiq. The 
upper bounds, on the other hand, can be calculated to order l/rriQ without additional HQET 
parameters. 

III. HADRONIC SIDE 

The hadronic matrix elements for semi-leptonic decay of a i? meson into a pseudoscalar 
meson P or a vector meson V may be parameterized as 

(V(v' e) \ Af^ \ B(v)) 

v' is the velocity of the final state meson, and the variable u = v ■ v' is a, measure of the 
recoil. One may relate uj to the momentum transfer g^ by cj = (M^ + M'jj — g^)/(2MBM/^). 
Therefore, with a proper choices of the current J^ and the four vector a^, one may readily 
select the form factor of interest and the corresponding sum rule for bounding. In the heavy 
quark symmetry limit, these form factors satisfy relations |jl|,0] 



h+{uj) = hviuj) = hA^{uj) = hA-,i^) = C{uj), 

h+{uj) = HaA^) = 0. (8) 

To bound Ha^ (or hy), one may choose J^ = A^ (or V^) and a^ = (0, 1, 0, 0). Then the 
factor to be bounded is "^^J |/iai(^)| (^^i5^ I^v(^)| ) and the sum rule used to bound 
is TiopE = TiHadronic- The Corresponding first excited state more massive than D* that 
contributes to the sum rule is the J^ = 1+ state, i.e. Di, since scalars do not contribute to 
Ti. 

For h+, one may take J^ = K^ and a^" = (1 + Ed/Md, 0, 0, -qs/Mo). Then the factor to 
be bounded is "^^ ' \h+{Lu)\ , and the combination of structure functions used in bounds is 

T(e) = -2 (1 + u) Ti + (1 + uf T2 + {Mb - Md - ef (1 + uf T4 

+2{MB-MD-e){l+LjfT,. (9) 

Since a'^ = v'^ + v'^, the first excited resonance that can contribute in this case is Di, due 
to the g^^"'^ structure of the D* form factor. 

Similarly, a convenient choice to isolate /i_ is to choose J^ = V^ and a^ = (1 — 
ED/MD,0,0,q3/Mn). Thus, 

T(e) = -2 (1 - cu)Ti + (1 - u;fT2 + {Mb + Md - tf (1 - uof T^ 

+2 {Mb + MD-t){l- ujf T5 (10) 

bounds '■'^~ -* |/i_(a;)| . -Di would be the first excited resonance for the same reason as in 
the case of /i+. 

It is impossible to single out /i^j and Ha-j, individually by any choice of a^. One good 
choice is to take J^ = A^ and a^ = {En/Mo, 0, 0, -Qs/Md). Then 

T(e) = -Ti + 00^X2 + {Mbu -M*j,-efT^ 

+2uj {Mbcu -M*B)-e) T5 (11) 

is the combination for bounding ^^^^^ |u;/ia2(^) + ^^3(^)1 • The first excited resonance that 
would contribute in this case would be the unobserved Dq. 



IV. PARTONIC SIDE 

Due to the heavy quark masses appearing in the problem, both the strong couphng 
constant asimg) (~ 0.3 at 2 GeV) and Aqqd /niQ will be good expansion parameters. Since 
uj is never very far from one, expanding in cj — 1 is also possible. We will keep terms up to 
order as{w — 1), dropping terms of order as{w — 1)^, a^, Aqcd/^^q ^^^ c^s-^qcd/itT'Q- 

The structure functions can be decomposed as 

rpfull ^j.l/m^j.a,^ (12) 

where the term Tj contains the tree and Aqcd/'^q contribution to the structure function, 
which has previously been calculated HoHTII. For these nonperturbative corrections, we keep 



the full u dependence. 

The leading order a^ corrections consist of bremsstrahlung radiation of a gluon from the 
heavy quarks and one loop virtual corrections. We expand to first order in w — 1, defining 

Tr = ^[U, + {w-l)V^. (13) 

TV 

The final results for the functions Ui and Vi are presented in the Appendix. The perturbative 
corrections to Ti were previously calculated in |T2|,p[], and agree with the results found here. 
We define the moments of the structure functions as 

^ fdee-T, = lt^ + A^\ (14) 

where 
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Ar = ^ fdee^^ [U, + {w-l) VI (15) 

ZTTI J IT 

The moments of Tj can be found in @,0. The moments for the perturbative corrections 
can be straightforwardly obtained from the functions in the Appendix. One thing to be noted 
is that the integration variable e in the bounds was defined in terms of hadronic variables. So 



when relating to the partonic computations, the corresponding integration variable should 
be changed to the one defined by partonic variables, namely, e^ = mi, — nic — v ■ q. The 
relation between them is eh = €p + 6, with 6 = Ec — EH + MB — rrib and Ec being the energy of 
the c-quark. We can now use these definitions to calculate the bounds on the form factors. 

V. BOUNDS ON INDIVIDUAL FORM FACTORS 

To form the bounds, one just takes the proper moments of the structure functions to 
form the combination required in the sum rules given in Section III. Corrections of order 
^qcd/"^Q5 af, as Aqcd/'^q, and a^ (tu — 1)^ should be small and have been neglected. To 
this order, the Aqcd/'^q corrections will depend on 3 HQET parameters; A, Ai and A2. From 
the measured mass difference, M/5. — M^, a very accurate value of A2 can be determined, 
A2 = 0.12 GeV^. It is very difficult to measure the parameters A and Ai individually, 
but a certain linear combination is much better determined |1^ . We will use three different 



parameter sets to show the dependence on A and Ai: (A) A = 0.4 GeV and Ai = —0.2 GeV^, 
(B) A = 0.3 GeV and Ai = -0.1 GeV^ (C) A = 0.5 GeV and Ai = -0.3 GeV^ 

The sum rule for bounding {u + 1)^ |/i^^(u;)| /(4c<j) uses Ti with the axial-axial current. 
The upper bound is simply the zeroth moment of Ti, which is by Eq. (|^) 



(^ + 1) 



2 



,(OMA 



/iAiHr</r""+Ar"". (16) 



4a; 

The variable A used for regularizing infrared divergences in the kinematic region away from 
zero recoil disappears in the final result. The first moment of Ti is needed for the lower 
bound, which is 

^^^ Ih^Mf > /f ^^^ + Af^^"" - ^ \ f/f ^^^ + A?^^) (17) 

4a; Ef)-^ — E^)* ^ ' 

The upper and lower bounds are shown in Fig. 1.0 For this section, the dotted curves are 



^For all the figures in this section we take m,\, = 4.8 GeV, rric = 1.4 GeV, a^ = 0.3 (corresponding 
to a scale of about 2 GeV), and A2 = 0.12 GeV^. The values of A and Ai are discussed in the text. 
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FIG. 1. Upper and lower bounds on {to + 1)^ |/i/i^(a;)| /(4c<j). The thick solid (dotted) curves 
are the upper and lower bounds including perturbative corrections for HQET parameter set (A) 
described in the text, and A = 1 GeV (2 GeV). The dashed curves are the bounds without pertur- 
bative corrections, also for HQET parameter set (A). The thin solid curves show the dependence 
on A and Ai, using parameter sets (B) and (C), with A = 1 GeV. 

the bounds without perturbative corrections using set (A) above, while the solid and dashed 
curves represent the bounds including the perturbative corrections with A = 1 GeV and 
A = 2 GeV, respectively. We have shown the bounds in the kinematic range 1 < cj ^ 1.25, 
where the higher order correction as{oo — 1)^ should be negligible. The thin solid curves use 
the other HQET parameter sets (B) and (C). 

A similar set of bounds for (cj^ — 1) 1/1^/(^)1 /(4a;) can be obtained readily by changing 
AA (axial-axial currents) to VV (vector-vector currents) in the above formulae. The bounds 
in this case are shown in Fig. 2. The tree level bounds go to zero at zero recoil because of 
the uj — 1 factor. 

The bounds on the other form factors involve higher moments of T4 and T5. The upper 



bound for {uo + 1)^ |/i+(cj)| /uo < ji^^"^ ^ where 
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FIG. 2. Upper and lower bounds on (o;^ — 1) |/iv'(u;)| /(4u;). The thick solid (dotted) curves are 
the bounds with perturbative corrections for HQET parameter set (A), and A = 1 GeV (2 GeV). 
The dashed curves are the bounds without perturbative corrections, also for HQET parameter set 
(A). The thin solid curves show the dependence on A and Ai, using parameter sets (B) and (C), 
with A = 1 GeV. 



7 upper 



■ UJ] 



1 



00 



-n 



(0)VV 



AO)VV j{2)VV 
-'2 "T -'4 



2{Mb - Md)II 



{l)VV 



\2 t{0)VV 



{Mb - MDyir + '^{Mb - MD)h 



{Q)VV 



2/. 



WW 



(18) 



l + uo 



Al 



(i))VV 



^{0)VV ^ j^(2)VV 



2{Mb - Md)A), 



{1)VV 



+ {Mb - M^)2Af ^^ + 2{Mb - Mr,)Af^^ - 2A«^^} . 



The lower bound is {uo + 1)^ |/i_|_(c<;)| /uo > W°^'^'^ , where 
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FIG. 3. Upper and lower bounds on (aj + 1)^ |/i-|-(ll')| /uj The thick solid (dotted) curves are the 
bounds with perturbative corrections for HQET parameter set (A), and A = 1 GeV (2 GeV). The 
dashed curves are the bounds without perturbative corrections, also for HQET parameter set (A). 
The thin solid curves show the dependence on A and Xi, using parameter sets (B) and (C), with 
A = 1 GeV. 



+ {Ms - MnfA^l^^"" + 2{Ms - Mz,)A«^'' - 2Af ^''l • 



They are plotted in Fig. 3. The perturbative physics pushes the bounds up from the tree 
level bounds in the region near the maximal momentum transfer while drags them down at 
large recoil. Changing A from 1 GeV to 2 GeV only slightly loosens both bounds. 

Without explicitly writing out the bounds on {uj — lY\h_{uj)\^ /uj, we simply state that 
they can be obtained from Eq. (|18|) and (|1^) by: 1) replacing each (1 +uj) factor by (1 —uj) 
and 2) substituting Mb — Mo by Mb + Mo. The bounds on this form factor are shown in 
Fig. 4. Notice that in this case, both the upper and lower bounds at tree level are identically 
zero. The perturbative corrections push both bounds away from zero. We still plot these 
negative lower bounding curves for reference even though the real lower bounds should be 
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FIG. 4. Upper and lower bounds on (a; — l)^|/i__(a;)p/u;. The thick solid (dotted) curves are the 
bounds with perturbative corrections for HQET parameter set (A), and A = 1 GeV (2 GeV). The 
dashed curves are the bounds without perturbative corrections, also for HQET parameter set (A). 
The thin solid curves show the dependence on A and Ai, using parameter sets (B) and (C), with 
A = 1 GeV. 

zero since the factors we are bounding are all positive-definite. Using A = 2 GeV in the 
calculation widens the upper bound by more than a factor of 2. 

Similarly, the upper bound for (u;^ — 1) |c<j/i^2(cj) + hA.j,(uj)\ /(4u;) < /iT^'^'", where 
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while the lower bound (cu^ — 1) |u;/i^2(u;) + hA.i{(jj)\ /(4tc;) < /i^^^'', where 
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FIG. 5. Upper and lower bounds on (w^ — l)\ujhA2{^) + hA.^{uj)\ /(4a-'). The thick solid 
(dotted) curves are the bounds with perturbative corrections for HQET parameter set (A), and 
A = 1 GeV (2 GeV). The dashed curves are the bounds without perturbative corrections, also for 
HQET parameter set (A). The thin solid curves show the dependence on A and Ai, using parameter 
sets (B) and (C), with A = 1 GeV. 
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Fig. 5 shows these bounds. The perturbative corrections drive the lower bound negative. 
Again, these lower bounding curves are plotted for reference purposes only. (The lower 
bound with A = 2 GeV is below the range plotted in this graph.) 

At (9(Aqcd/'"^q), the upper bounds will not depend upon Ai and A2. However, since 
l/(M£)i — Mb) ~ 1/Aqcd they will affect the lower bounds. It is possible to obtain the 
upper bounds to order 0{K^Qj^/m^Q), since at this order there are no new parameters in the 
OPE. These corrections only slightly modify the above upper bounds, so we do not show 



13 



them here. 

VI. COMPARISON WITH MODELS 

We choose from the hterature the following commonly used form factor models for com- 
parison with our bounds: 



;i) ISGW2 Model ]14|,|15 



(2) Light-Front Model ^^ 



(3) BSW II Model [|lf^ 



(4) NS Model |21 



(5) COQM Model Q. 

The ISGW2 model is the updated version of the original ISGW model which incorporates 
the constraints of heavy quark symmetry and relativistic corrections. The BSW II model, 
derived from the BSW I model by improving the pole structure of form factors, is considered 
with updated pole masses ^^. The NS model was proposed as a simplified alternative to 
the NRSX model [^ which specifies the Isgur-Wise function by several strong assumptions 
such as the pole structure and the condition for the derivative of the single-pole form factor. 
The form factors are then obtained by employing the heavy quark symmetry relations. The 
COQM model uses a covariant oscillator quark model to calculate the Isgur-Wise function, 
and then uses heavy quark symmetry to relate it to the form factors. 

Figs. 6-10 are plots for the different form factors from the models and the bounding 
curves. In plotting these figures, we used nib = 4.8 GeV, rric = 1.4 GeV, as = 0.3 (corre- 
sponding to a scale of about 2 GeV), A2 = 0.12 GeV^, A = 1 GeV and current PDG data 
on heavy meson masses. As mentioned above, while A2 can be measured from the mass 
difference M^* — M^ to a high accuracy, Ai and A are not easy to obtain experimentally. 
Here we pick the parameter set (A) discussed above, A = 0.4 GeV, Ai = —0.2 GeV^. This 
uncertainty on Ai and A will slightly modify our bounds to the order and kinematic range 
we are working, as seen from Figs. 1-5. 

14 



0.9 



3 



0.8 



< 



^0.7 
3 



0.6 



0.5 




1.25 



FIG. 6. The model values of {lo + 1)^ |/i^j(a;)| /(4a;) along with the corresponding bounds for 
comparison. The thick solid lines are the upper and lower bounds. The thin solid curve is the 
ISGW2 model. The long dashed curve is the LF model, and the short dashed curve is the BSW II 
model. The dot-dashed curve is the NS model. The dotted curve is the COQM model. 

Fig. 6 shows the model values of {u + 1)^ |/iyi^(cj)| /(4cj) along with the corresponding 
bounds for comparison. The first three models have curves falling within or close to the 
perturbative bounds, while the NS and COQM models have curves slightly off the bounds 
near zero recoil because they are models designed for the Isgur-Wise function without taking 
into account the perturbative corrections. The COQM model, however, has a curve which 
falls far below the lower bound at large ui. 

In Fig. 7, the ISGW2 and NS models lie above the upper bound for (a;^ — 
1) \uhji^ + Ha^I /(4cj) at large recoil, while the other models are closer to our bounds. 
For A = 1 GeV they are above our bounds for most of the kinematic range. If we take 
A = 2 GeV, however, they would be within our bounds. (ISGW2 and NS would still be too 
large.) 
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FIG. 7. The model values of (o;^ — 1) \ujhA2 + ^^3! /(4u;) along with the corresponding bounds 
for comparison. The curves are labeled the same as in Fig. 6. 

Only the 1SGW2 model predicts a curve near our bounds for (cj^ — 1) \hv\ , shown in 
Fig. 8; the rest are too small. The Light-Front and BSW 11 models will be between our 
bounds when the scale is set at A = 2 GeV , but the NS and COQM models still fall below 
the lower bound. 

As shown in Fig. 9, the 1SGW2 model agrees with our bounds for {u + 1)^ \h^\ juj very 
well. For the same reason as in Fig. 6, the NS and COQM models start from the position 
predicted by heavy quark symmetry at zero recoil, without perturbative corrections. The 
Light-Front, NS, and COQM models stay below the lower bound, and the BSW 11 model 
starts above the bounds at c<j = 1, then cross the bounds and enter the region under the 
lower bound. At large recoil, the 1SGW2, BSW 11, and NS model predictions would lie 
within the bounds for most of the kinematic range with a larger value for A. 
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FIG. 8. The model values of {u"^ — 1) \hy\ along with the corresponding bounds for comparison. 
The curves are labeled the same as in Fig. 6. 

Fig. 10 shows the curves for {uo — 1)^ |/i_| /cj, where most models are consistent with 
our bounds within the kinematic range of interest. The curve for the COQM model is zero, 
since it uses heavy quark symmetry, which gives h_{w) = in the infinite mass limit. 

The scale we choose to plot these diagrams is A = 1 GeV, since this gives the tightest 
bounds. Had we chosen 2 GeV as our working scale, the bounds would be much less stringent 
and thus would accommodate the models which originally fell slightly outside our bounds 
in these diagrams. 

VII. ORDER a|/3o CORRECTIONS AT ZERO RECOIL 



By using the techniques introduced in Ref. |^, we can calculate the a^Po contribution 
to the structure functions. Typically, these corrections are about 90% of the full a^ rate, so 
they can be used as a rough estimate of the next order corrections. The O^a^Po) correction 
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FIG. 9. The model values of {uj + 1)^ \h+\ /^ along with the corresponding hounds for compar- 
ison. The curves are labeled the same as in Fig. 6. 

to the structure functions, T"^ , can be related to the 0{as) correction, T^"{X)^ calculated 
with a gluon mass A, as [E^] 



r^^° . _^^^lli^ r ^ f^«.(A) - ,^7;-(0)^ 



(22) 



An Jo X'^ \ ' ' ' A2 + A2 * / 
where a^^^ is the strong coupling constant evaluated in the "V-scheme" [^, and is related 
to the MS coupling constant a^ by 



a7\f^)=Mf^) + m^Po + 



(23) 



3 47r 

We numerically calculate these corrections at zero recoil. In analogy to Eq. (|I4|) , we write 
the corrections to the moments of the structure as 



— /rfee"T, = 4"^ + A^") 
2m 



(24) 



where Af" now contains both the 0{as) and the 0(a^/3o) contribution to the n^^ moment 



of the i structure function at zero recoil. The results for the A) s for the zeroth and first 
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FIG. 10. The model values of (w — 1)^ |/i_| /u> along with the corresponding bounds for com- 
parison. The curves are labeled the same as in Fig. 6. 

moments are presented in Table 1, for A = 1,2 GeV. The 0{a1l3o) contributions are often 
large, especially for A = 1 GeV. This may be an indication that we need to go to larger A 
to get a reliable perturbative expansion. 

If we look at the actual combination of structure functions that appear in our bounds, 
the situation is more promising. The 0(a^/9o) corrections to the upper and lower bounds 
on the form factors at zero recoil are shown in Table 2 for A = 1 GeV. In this case, the 
0(a^/3o) corrections tend to be rather small, so it seems that the perturbative expansion for 
the bounds is under control. 

VIII. DISCUSSION 

As discussed earlier, our bounds on the form factors are derived to order Aqcd/'"^q and 
contain leading perturbative QCD corrections. They are constructed to be model indepen- 
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dent and can be used to test various model predictions for the form factors. Furthermore, 
they may be used to bound inclusive or exclusive decay rates. The validity of the lower 
bounds requires the extra assumption of negligible multiparticle production in the decays. 
This assumption is rather mild, since it is true in the large Nc limit and is supported by 
current experimental data. The upper bound is rigorously valid without this extra assump- 
tion. Therefore, any significant deviation of the phenomenological form factor beyond the 
upper bound indicates the need for some modification of the model. In general, our bounds 
are stricter and more accurate near zero recoil. 

The bounds become much less stringent when A is increased. Therefore, we should use 
the smallest value of A for which there is a reasonable perturbative expansion. We also 
observed that the 0(q;^/9o) corrections to the bounds at zero recoil are small for A = 1 GeV, 
which in turn suggests that a perturbative expansion in this problem seems to work well, 
provided that the 0(ag/9o) corrections dominate in the complete 0{a1) corrections. 

Among the models we considered in this work, the ISGW2, LF, and BSW II models 
"pass" most of the qualifications with A = 2 GeV. For smaller values of A, ISGW2 seems 
to do the best. However, all models fail, over essentially the whole kinematic range, for 
the form factor combination ^^^-^ \0ohA2 + ^^3! , and thus the models should be modified to 
satisfy this bound. 
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APPENDIX: RADIATIVE CORRECTIONS TO THE STRUCTURE FUNCTIONS 

The structure functions from perturbative corrections are written in the form Tf'^ = 
{as/'n')[Uf'^ + {uj — 1) V/'^] for the ith structure function and current J. rrib is the mass of 
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the 6-quark, rric is the mass of the c-quark, z = mc/mb, u = v ■ v', and A is the mass of the 
gluon used for infrared regularization. We have taken A ^ wherever possible. 
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A = 2 GeV 
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TABLE I. Order Og and ajj3o corrections to the moments of the structure functions at zero 
recoil, evaluated for A = 1, 2 GeV. For this table, as = 0.3, and Uf = 4. 
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Tree Level 0{as) ©(a^^o) 


Upper (^|/i^j2 
Lower (^|/^^J2 


1 -0.071 0.016 
1 -0.084 0.006 


Vpper^\u;hA,+hAf 
Lower ^\u;hA^+hAf 


0.011 0.007 
-0.005 -0.001 


Upper ^\hv\^ 
Lower ^|Vp 


0.006 0.004 
-0.003 -0.001 


Upper (^|/.+ |2 
Lower ("+^)V+P 


4 0.20 -0.05 
4 0.19 -0.06 


Upper ("-^)>_p 
Lower ("-/)' |/i_p 







TABLE IL Bound on form factors at zero recoil, evaluated with A = 1 GeV. 
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